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Abstract 

A crucial point in statistical mechanics is the dehnition of the notion of thermal equilibrium, which 
is given as the state that maximises the von Neumann entropy, under the validity of some constraints. 
Arguing that such a notion can never be experimentally probed, in this paper we propose a new notion 
of thermal equilibrium, focused on observables rather than on the full state of the quantum system. We 
characterise such notion of thermal equilibrium for an arbitrary observable via the maximisation of its 
Shannon entropy and we bring to light the thermal properties that such a principle heralds. The relation 
with Gibbs ensembles is studied and understood. We apply such a notion of equilibrium to a closed 
quantum system and show that there is always a class of observables which exhibits thermal equilibrium 
properties and we give a recipe to explicitly construct them. We bring to light an intimate connection of 
such a principle with the Eigenstate Thermalisation Hypothesis and provide the first prediction about 
which observables satisfy it. 


Introduction 

To understand under which conditions thermodynamics emerges from the microscopic dynamics is the ulti¬ 
mate goal of statistical mechanics. However, despite the fact that the theory is more than 100 years old, we 
are still discussing its foundations and its regime of applicability. The ordinary way in which thermal equilib¬ 
rium properties are obtained, in statistical mechanics, is through a complete characterisation of the thermal 
form of the state of the system. One way of deriving such form is by using Jaynes prOTcjpZe[T][21[31ll], which 
is the constrained maximisation of von Neumann entropy S'vn = Trplogp. The main point of Jaynes work 
was that statistical mechanics can be seen under the light of probabilistic inference, in which one is forced 
to give prediction about some macroscopic properties of the system, even though the information that we 
have on the system is not complete. Adopting such point of view, Jaynes showed that the unique state that 
maximises 5 'vn (compatibly with the prior information that we have on the system) is our best guess about 
the state of the system at the equilibrium. The outcomes of such procedure are the so-called Gibbs ensembles. 

In the following we argue that such a notion of thermal equilibrium, de facto is not experimentally 
testable because it gives predictions about all possible observables of the system, even the ones which we are 
not able to measure. To overcome this issue, we propose a weaker notion of thermal equilibrium, specific for 
a given observable and therefore experimentally verifiable. 

The issue is particularly relevant for the so-called “Pure states statistical mechanics” 0 0 [3 m [a [13 
nmnnKHnnKiKniiiKis], which aims to understand how and in which sense thermal equilibrium 
properties emerge in a closed quantum system, under the assumption that the dynamic is unitary. In the 
last fifteen years we witnessed a revival of interest in these questions, mainly due to remarkable progresses 
in the experimental investigation of isolated quantum systemsIlHl [21] [22 [23l [24] [25] ■ Nowadays, the high 
degree of manipulability and isolation from the environment that we are able to reach makes possible to 
experimentally investigate such kind of questions and to probe the theoretical predictions. 


The starting point of Jaynes’ derivation of statistical mechanics is that i^vn is a way of estimating the 
uncertainty that we have about which pure state the system inhabits. Unfortunately, even though 5 'vn is 
undoubtedly an important quantity, we know from quantum information theory that it does not address all 
kind of ignorance we have about the system. In, it is not the entropy of an observable (though the state is 
observable); its conceptual meaning is not tied to something that we can measure. 

This issue is intimately related with the way we acquire information about a system, i.e. via measure¬ 
ments. The process of measuring an observable O on a quantum system allows to probe only the diagonal 
part of the density matrix (A^lplAi), when this is written in the observable eigenbasis {|Ai)}. For such a 
reason, from the experimental point of view, it is not possible to assess whether a many-body quantum 
system is at thermal equilibrium (e.g. Gibbs state pc)- the number of observables needed to probe all the 
density matrix elements is too big. In any experimentally reasonable situation we have access only to a few 
(sometimes just one or two) observables. It is therefore natural to imagine situations in which the outcomes 
of measurements are compatible with the assumption of thermal equilibrium, while the rest of the density 
matrix of the system is not. 

Despite that, we think that the fact that a distribution is compatible with its thermal counterpart will 
lead to the emergence of certain thermal properties, concerning the specific observable under scrutiny. Build¬ 
ing our intuition on that, we propose a new notion of thermal equilibrium specific for a given observable, 
experimentally verifiable and which relies on a figure of merit that is not the von Neumann entropy. A good 
choice for such a figure of merit comes from quantum information theory and it is the Shannon entropy 
Ho of the eigenvalues probability distribution {p(Aj)} of an observable O. The well-known operational 
interpretation of Ho^^ matches our needs since it addresses the issue of the knowledge of an observable 
and it provides a measure for the entropy of its probability distribution. 

Throughout the paper we will work under the assumption that the Hilbert space of the system has finite 
dimension and we will refer to the cases in which the Hamiltonian of the system has no local conserved 
quantities, even though it is possible to address situations where there are several conserved quantities, like 
integrable quantum systems. Moreover, we will assume that the observable has a pure-point spectrum with 
the following spectral decomposition O = where Hj is the projector onto the eigenspace defined 

by the eigenvalue Xj. Hq is the entropy of its eigenvalues probability distribution p{Xj) = Tr {pHj) 


Ho[p] = - 

j ■■ AjGcrei 
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where ao is the spectrum of O. 

We propose to define the notion of thermal equilibrium, for an arbitrary but fixed observable O, via a 
characterisation of the probability distribution of its eigenvalues. We will say that 

O is at thermal equilibrium when its eigenvalues probability distribution p{Xj) maximises the Shannon 
entropy Hq , under arbitrary perturbations with conserved energy. We call an observable with such a proba¬ 
bility distribution: thermal observable. 

It is important to note that such a principle characterises only the probability distribution at equilibrium 
and it does not uniquely identify an equilibrium state. Given the equilibrium distribution p{Xj) = p®‘^(Aj), 
there will be several quantum states p which give the same probability distribution for the eigenvalues Xj. 
In this sense this is weaker a notion of equilibrium, with respect to the ordinary one in which the whole state 
of the system is required to be an equilibrium state. 

In the rest of the paper we study the main consequences of the proposed notion of observable-thermal- 
equilbrium: its physical meaning and the relation with Gibbs ensembles. These results provide evidences that 
the proposed notion of equilibrium is able to address the emergence of thermalisation. This is our first result. 
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Furthermore, we study the proposed notion of equilibrium in a closed quantum system and prove that 
there is a large class of bases of the Hilbert space which always exhibit thermal behaviour and we give an 
algorithm to explicitly construct them. We dub them Hamiltonian Unbiased Bases (HUBs) and, accordingly, 
we call an observable which is diagonal in one of these bases Hamiltonian Unbiased Observable (HUO). The 
existence and precise characterisation of observables which always thermalise in a closed quantum system is 
our second result. Furthermore, we investigate the relation between the notion of thermal observable and 
one of the main paradigms of pure states statistical mechanics: the Eigenstate Thermalisation Hypothesis 
(ETHipTl dilMlinnilSIlISlI^IMlISSlIMlISZlIMlIMlIlS]- We find an intimate connection between the 
concept of HUOs and ETH: the reason why these observables thermalise is precisely because they satisfy the 
ETH. Hence, with the existence and characterisation of the HUOs we are providing a genuine new prediction 
about which observables satisfies ETH, for any given Hamiltonian. The existence of this relation between 
HUOs and ETH is a highly non trivial feature and the fact that we can use it to predict which observables 
will satisfy ETH is our third result. Indeed, as far as we know, there is no way to predict which observables 
should satisfy ETH in a generic Hamiltonian system. It is an hypothesis which has to be checked by inspec¬ 
tion, case by case. 

In the conclusive section we summarise the results and discuss their relevance for some open questions. 


Results 

Information-theoretic equilibrinm 

The request that the equilibrium distribution must be a maximum for Ho can be phrased as a constrained 
optimisation problem, which can be solved using the Lagrange multiplier technique. The details of the 
treatment can be found in the Methods section. Two sets of equilibrium equations are obtained and we show 
in which sense they account for the emergence of thermodynamic behaviour in the observable O. We assume 
that the only knowledge that we have on the system is the normalisation of the state and the mean value of 
the energy (T) = Eq, where T is the Hamiltonian of the system. We also make the following assumptions 
on the Hamiltonian: that it has the following spectral decomposition T = EaTa, where = \Ea) {Ea\ 
and that its eigenvectors {|Eq,)} provide a full basis of the Hilbert space. We call |i/'n) the eigenstates of the 
density operator, pn are the respective projectors and its eigenvalues. To describe the state of the system 
we use the following convenient basis: {|j, s)} in which the hrst index j runs over different eigenvalues Xj 
of O and the second index s accounts for the fact that there might be degeneracies and it labels the states 
within each subspace defined by A^. We also make use of the projectors H^s = \j,s) (j, s|. Eurthermore, 
we are going to need the following quantity £n{j,s) = (ipn\^j,sT Itpn), which has the following physical 
interpretation. Hj ^T is the projection of the Hamiltonian onto the eigenstate |j, s). Hence we understand 
£n{j,s) as the “average (on \tpn)) amount of energy stored in the basis element |j, s)”. Using the overlaps 
= (j, sIV’ra) the expression of the constraints: 


Cn = Tr(p) - 1 = gn 



2 

- 1 


n\j,s 


Ce = Tr {pT) - Eo = y] qn£n{j, s) - Eq . 

n-,j,s 


(2a) 

(2b) 


It is important to remark here that the knowledge of Eq is always subject to uncertainty, which we call 
6E. In any physically reasonable situation there will be two conceptually different sources of uncertainty: 
a purely experimental one and an exquisitely quantum one. In this sense, all the states p : TrpT € Iq = 
[Eq — ^,Eo + will be solutions of the constraint equation, up to a certain uncertainty SE. Even though 
we do not make any assumption on such a quantity, we note that 6E is usually assumed to be small on a 
macroscopic scale but still big enough to host a large number of eigenvalues of the Hamiltonian. 
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Exploiting Lagrange’s multipliers technique we obtain four sets of equations. Derivatives with respect to 
the multipliers enforce the validity of the constraints while the derivatives with respect to the field variables 
give two independent set of equations. Using two particular linear combinations of them we obtain the 
following equilibrium equations (EEs): 


s) — s) 


I E P I = (1 - - XE£nU,s) , 


(3a) 

(3b) 


where Xe and Xn are the Lagrange multipliers associated to C^; and Cn, respectively. 

Our first result is the understanding of the physical consequences of these equations. Together, they 
account for the emergence of thermal behaviour in O. The first one l (l3all l gives the stability under the flow 
generated by the Hamiltonian and it implies that the equilibrium distribution peq{Xj) has to be be invariant 
under the unitary dynamics. Indeed, writing the time evolution equation for p{Xj), using the von Neumann 
equation, we obtain 


*S^p(Aj) = E {£n{j, s) - £n{j, s)) = 0 , (4) 

n,s 

where the superscript “eg” stands for “at the equilibrium”, i.e. after plugging in the EEs. 

The second equation l (l3bll l fixes the functional form of the distribution with respect to the Hamiltonian 
and to the Lagrange multipliers. It can be shown that it is responsible for the emergence of a thermody¬ 
namical relation between Ho and the mean value of the energy. Integrating (I3bll over the whole spectrum 
of O we obtain 


= {1 - Xn) - Xe Eo. (5) 

There is a linear contribution in the mean value of the energy, plus a “zero-point” term = (1 — 
Xn)- This relation brings to light the thermodynamical meaning of Shannon entropy since such linear 
dependence on the average energy is a distinguishing feature of thermodynamic equilibrium. We note a 
strong analogy with the properties of von Neumann entropy, which acquires thermodynamical relevance 
once the state of the system is the Gibbs state pc = ^ ’ 'where Z = Tre“^^ is the partition function 


SMpg) = log Z + pEo ^ = (1 - Ajv) - Ab Uo . 


( 6 ) 


Relation with statistical mechanics 

We now come to the issue of the relation with the ordinary notion of thermal equilibrium and, therefore, 
with the Gibbs ensemble. With respect to ours, this is a much more stringent condition because it is a 
characterisation of the full state of the system. For such a reason, in order to understand if our notion of 
thermal observable is compatible with the standard characterisation of thermal equilibrium, we need to find 
the condition under which our criterion gives a complete characterisation of the state of the system. Con¬ 
sidering that we are using a maximum-entropy principle, a plausible auxiliary condition is the maximisation 
of the smallest among all the Shannon entropies. Such a request fully characterises the state of the system 
because the lowest Shannon entropy of the state is unique and it is the one in which the density matrix of 
the state is diagonal. Indeed, using the Schur-concavity of the Shannon entropy and the Schur-Horn theorem 
on eigenvalues of an Hermitian matrix |26j it is easy to prove that 


{eq:EEl} 

{eq:EE2} 


{eq:von} 
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min Ho{p) = S^n{p) = -Tr {plogp) , (7) 

O&A 

where A is the algebra of the observables [H]. Our minimalist request to maximise the lowest Shannon 
entropy is translated in the condition of maximising von Neumann entropy. The state which maximises von 
Neumann’s entropy, with the constraint of fixed total energy, is precisely the Gibbs state. It is therefore clear 
that our proposal constitutes an observable-wise generalisation of the ordinary notion of thermal equilibrium. 

In the next section we will apply the proposed notion of thermal observable to the so-called “pure-states 
statistical mechanics” [3 El El 0 0 llOj . which aims at understanding when, and in which sense, closed 
quantum system might come to thermal equilibrium. In particular we will investigate the relation with the 
main paradigm of pure-states statistical mechanics: the ETH. 

Closed Quantum Systems - Relation to ETH 

ETH, in its original formulation [MUSniEIlEaiMllMlESlEe], is an ansatz on the matrix elements of an 
observable when it is written in the Hamiltonian eigenbasis {jifa)}: 


^oETH 

^a/3 


f^o\E)So.p + e-^fg\E,oj)R^f, , 


( 8 ) 


where E = uj = E^ — Ep while and are smooth functions of their arguments. S{E) 

is the thermodynamic entropy at energy E defined as = E'^^Se{E — Ea)., where (5^ is a smeared 

version of the Dirac delta distribution. Rap is a complex random variable with zero mean and unit variance. 
Furthermore, it is important to remember that ETH by itself does not guarantee thermalisation, we need 
to impose that the initial state has a small dispersion in the energy eigenbasis [35] • When this is true one 
says that thermalises in the sense that its dynamically evolving expectation value is close to the 

microcanonical expectation value defined by the conditions on the average energy 


(9) 

where pmc = PmciEo, 5E) is the microcanonical state defined by the condition on the average value of the 
energy TrpT G Iq. 

The reason why Eq. ([9]) becomes true if ETH holds is purely technical, in the following sense. The small 
dispersion assumption on the initial state is telling us that when we expand the initial state in the energy 
eigenbasis, we will only have non-zero contributions from a certain energy window: 


\i>o) ='^Ca\Ea) ^ ^ Ca\Ea) ■ ( 10 ) 

a'.Eoi^Io 

If such dispersion SE is small enough, the smooth functions \ /q^ will not vary too much in such an 
energy window. Once ETH has been checked by inspection, one can conclude that there is thermalisation 
of the average value in the sense of Eq.Q. It is important to note that, from the conceptual point of view, 
the “small energy-dispersion assumption” is a key element of the emergence of thermal equilibrium but it 
has nothing to do with ETH which, by itself, is only the ansatz in Eq.®. Nevertheless, this assumption is 
expected to hold in real experiments because, when working with a many-body quantum system, it is almost 
impossible to prepare coherent superpositions of states with macroscopically different energies [HI EH. 

The main intuition behind the emergence of ETH for an observable O (“Berry conjecture” [53EH1) is 
that for a many-body quantum system the Hamiltonian eigenstates should be so complicated that, when 
they are written in the eigenbasis of O, their coefficients can be effectively described in term of randomly 
chosen coefficients and this would lead to the emergence of thermalisation. So far, nothing has been said 


{eq:ETH} 


{eq:AverTherm} 


{eq:pure} 
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about which observables we should look at, however it is generally expected that such a statement should 
be true for local observables. We explicitly show that such mechanism is revealed by the maximisation of 
Shannon entropy and this allows us to give a prediction about which observables should satisfy ETH. In the 
next section we will also clarify why we expect local observables to satisfy ETH and we will also provide 
some general conditions on the Hamiltonian whose validity guarantees that local observables will satisfy ETH. 

Before we continue we need to define the following short-hand notation: ct G Iq 

stands for a : G Iq- 

Hamiltonian Unbiased Observables and ETH 

To study the relation with ETH we need to change our perspective. The point of view that we are adopting 
in the following. One of the key-points behind the ETH is that, in many real cases, the expectation values 
computed onto the Hamiltonian eigenvectors can be very close to the thermal expectation values. Moreover, 
when one wants to argue that thermalisation in a closed quantum system arises because of ETH, a main 
assumption is that the initial pure state of the system I'i/jg) bas a very small energy uncertainty AE, with 
respect to the average energy Eq: <C 1. Following these two insights, we take the extreme limit in which 

AE = 0. With such a choice, we are left with an Hamiltonian eigenstate and a constraint equation given 
by TipH = En- We note that Ea. (l5a1) is trivially satisfied for an Hamiltonian eigenstate. Hence we assume 
that |'0) = \Ea) : Ea G Iq and use the solvability of Ea. (l3bl) as a criterion to look for observables which can 
be thermal. While this is a very specific choice, we will show that it is enough to unravel some interesting 
feature regarding the ETH. 

With such an assumption, the second equilibrium equation becomes 


- log = (1 - Av) - Xe\D,\^E^ (11) 

Dividing by \Dj\'^ on each side we obtain a right-hand side which does not depend on the label j. Hence, 
keeping in mind that xlogx -G 0 for x —> 0, the most general solution of this equation is a constant 
distribution with support on some subset (la C oq depending on Ea) of the spectrum and zero on its 
complementary: 

[ ^ ^ XjGXaCao 

P{^j) ~ ^ (12) {eq:const} 

[ 0 V Aj G OolXa 

where da = dimla is the number of orthogonal states on which the distribution has non-zero value. The 
distribution of eigenvalues of O given by Eg. (1121) fully agrees with the prediction from the microcanonical 
ensemble Pmc, defined by the condition in Eg. (1121) . This is true, in particular, for the expectation value 


(O) 


■E 

j :Xj£Xa 


A, = Tr 




(13) 


It has to be understood here that Ea. dT^ is a highly non-trivial condition on O, which is not going to be 
fulfilled by every observable since it imposes a very specific relation between its eigenstates and the energy 
eigenvectors. 


Within the context of quantum information theory, and especially in quantum cryptography, a similar 
relation has already been studied. Two bases of a T>—dimensional Hilbert space (By = {|r)i)} and B^ = 
{|?Cj)}) are called mutually unbiased bases (MUB) [42l|43l|44] when 


= X V i,j = 1,...,T> (14) 

Such a concept is a generalisation, expressed in term of vector bases, of canonically conjugated operators. 
In other words, each vector of By is completely delocalised in the basis By, and viceversa. Here we mention 
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the result about MUBs which matter most for our purposes: given the 2^-dimensional Hilbert space of N 
qubits, there are 2-^ + 1 MUBs and we have an algorithm to explicitly find all of them[5^. Therefore, if 
we are in an energy eigenstate, an observable unbiased with respect to the Hamiltonian basis will always 
exhibit a microcanonical distribution. This is also true if our state is not exactly an energy eigenstate, it is 
enough to have a state that has a sufficiently narrowed energy distribution. We now provide a simple argu¬ 
ment to prove such a statement. We also note that such condition is closely related to the small dispersion 
condition briefly discussed before and that it is necessary to guarantee thermalisation, according to the ETH. 

By assumption, the pure state in Ea. dTOl) has an energy distribution {pa = with small disper¬ 

sion. This implies that its Shannon entropy Ht has a small value, because the profile of the distribution 
will be peaked around a certain value. For such a reason Ht{\iPo)) will be much smaller than its maximum 
value 


Ht{\'^o)) < logT> 


(15) 


Moreover, it can be proven that between any pair of MUB, like the Hamiltonian eigenbasis and an HUB, 
there exists the following entropic uncertainty relation, involving their Shannon entropies: 


(I</>))+I</>))> log!? V|0)sH 


(16) 


Putting together Eg. dTSl) and Ea. dTBl) we obtain that, for all the states with a small energy dispersion, 
the Shannon entropy of an HUB will always have a value close to the maximum: 


HHUB{\i^o))-\ogV (17) 

This in turn implies that the distribution of all the HUO will be approximately the same as the one 
computed on the microcanonical state. 

For such a reason, we now study the properties of a HUO: 


= I-?’ 


{j,s\Ea) 


Vv 


(18) 


HUOs and ETH: diagonal matrix elements 

In order to investigate the relation with ETH we need to study the matrix elements of a HUO in the energy 
basis: 


with tLi“/ = { 0 js ,/3 — Qjs,a)- It is straightforward to conclude that its diagonal matrix elements are 
constant in such a basis and therefore the so-called (Hamiltonian) Eigenstate Expectation values reproduce 
the microcanonical expectation values: 

Oaa = ^ E = ^Tr O = Tr (O p„,) = (20) 

3,s 

This is the first part of our third result and as it is, it can already be used to explain the emergence 
of thermalisation in closed quantum system. In m US] Reimann proved an important theorem about 
equilibration of closed quantum systems. He was able to show that under certain experimentally realistic 


{eq:SmallH} 


{eq:EntUnc} 


{eq:Const} 


{eq:off} 


{eq:EE} 
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conditions, the mean value of an observable is not much different from its value computed on the time- 
averaged density matrix, or Diagonal Ensemble (DE): 

PDE = ''^\Ca\^ \Ea) {Ea\ (21) {eq:DE} 

OL 

where Cq, = (V'oI^'q) and \^po) is the initial pure state of the isolated system. Roughly speaking, the two main 
assumptions made by Reimann are the following: first, that in any experimentally realistic condition the 
state of the system will occupy a huge number of energy eigenstates, even if the average energy is known up 
to a small experimental uncertainty; second, that the observable under study has a finite range of average 
values, due to the fact that we wish to measure it. For a clear and synthetic discussion on this topic we 
suggest [14] and we send the reader to the original references [mill]. We note that the first assumption 
does not contradict the small energy dispersion assumption. Indeed, as argued by Reimann, in a many-body 
quantum system, even if the energy is known up to a macroscopically small scale 6E, there will be a huge 
number of eigenstates within the range E^ G Iq- This is formally equivalent to say that we will observe a 
high density of energy eigenstates. To conclude, given a HUB it is always possible to obtain a HUO which 
satisfies the finite-range assumption. We can therefore apply Reimann’s theorem to HUOs. 

It is important to note that Reimann’s theorem explains equilibration around the DE but this does 
not necessarily entail thermalisation. The DE still retains information about the initial state while, on the 
contrary, thermalisation is defined (also) by the independence on the initial state. This is the point where 
our result is able to take a step forward and explain the emergence of thermal equilibrium in the HUOs. We 
can use Ea. (|20)l to prove that all the HUOs exhibit complete independence from the initial conditions: 


Tr (Opde) = Y. = Tr {Opmc) = 

O' 


( 22 ) 


HUOs and ETH: off-diagonal matrix elements 

In order to prove that a HUO satisfies ETH we need to study also its off-diagonal matrix elements. By 
using Eg. ([151) . we can investigate how the phases = (%s ,/3 — (^js,a) are distributed. This can be done 
numerically, exploiting the available algorithms to generate MUBs[^. The numerical investigation of the 
distribution of is reported in the supplementary material. Here we simply state the result: for a fixed 
value of the energy quantum numbers, the observed distributions of cosa;“/, sinw“j^ are well described by 

the assumption that are independent and randomly distributed in [— 7r,7r], with a constant probability 
distribution. 

There are different ways in which this result can be used. The general argument is the following and 
it reflects the spirit of ETH: in Ea. (fT9l) . the phases will have a randomising action on the eigenvalues 
Xj and this will make the value of the off-diagonal matrix elements severely smaller than the value of the 
diagonal ones: 


< ^^rO = Oaa (23) 

The randomness of the coefficients involved in the evaluation of the off-diagonal matrix elements has 
been recently proposed as the basic mechanism to explain the scaling behaviour which was observed to 
occur in some models [SU EH [531 [Ml 155] . The maximisation of Shannon entropy is therefore giving us a 
recipe to find the observables for which this is true. 


HUOs and ETH: two important examples 

In order to understand how this works in practice we need to say something specific about the eigenvalues. 
We highlight two important cases in which our result is helpful: an observable which is highly degenerate 
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and an observable whose eigenvalues distribution is not correlated with the phases . 

Highly degenerate observable - Assuming that j = 1,..., while s = 1,..., 1^2 with T? ^ ^ 'Hi, the 

sum in Eq. (HID splits into Hi of sums and each one of them is a sum of (D 2 ^ 1 identically distributed 
random variables. We can apply the central limit theorem to the real and imaginary part of Eq. dUD and 
obtain the following expression for the off-diagonal matrix elements 


ReOty 


■Di 

AjA/’j 

f=l 


VH2 J Vh 


TZn 


(24) 


In which indicates a gaussian probability distribution with mean x and variance cr^ and TZap 

is a zero mean and unit variance random variable. This proves that a highly degenerate observable satisfies 
the ETH ansatz, Eq.(|8D. The diagonal matrix elements reproduces the microcanonical expectation values 
and the off-diagonal matrix elements are well described by a random variable with zero mean and ^ 
variance. This is precisely what we expect from a local observable, since its eigenvalues have a huge number 
of degeneracies, which grows exponentially with the size of the system, and it is in full agreement with the 
randomness conjecture made in |541155) . 


Uncorrelated distribution - If the eigenvalues distributions and the phases | are not 

correlated Ea. dTOl) becomes 


Where we used the fact that the two sequences are uncorrelated to approximate the sum as the product 
of two sums and in the second identity we used the fact that {|j, s)} is a complete basis. We conclude that 
the off-diagonal matrix elements of such an observable are much smaller than the diagonal ones and therefore 
we can neglect them when we compute its dynamical expectation value: 

a,/3 a 

The assumption that the sequences {^j}j g | are not correlated will not be always fulfilled, 

nevertheless this is what we intuitively expect from the following reasoning. We know that a;“/ can be 
described by a random variable with a constant probability distribution, this suggests that there is no corre¬ 
lation between |a;“/1 and its labeling {j, s}^ g. On the contrary, the eigenvalues Xj are highly correlated 

to the labels and therefore we do not expect them to be correlated to the |w“/1 . 

These results follows from the study of the equilibrium equations, under the assumption that the state 
is an Hamiltonian eigenstate IV') = \Ea) belonging to the energy shell Iq- It is straightforward to see that 
the same results hold when the state is the microcanonical state Pmc{Eo, SE) involved in Eq.([^ and defined 
by the condition TrpT G Iq. 

In the next section we will exploit the proven results to investigate when local observables are expected 
to satisfy ETH, in a closed quantum system. 


Discussion 

We proposed a new notion of thermal equilibrium for an observable O: we say that O is a thermal observable 
when its eigenvalues probability distribution maximises its Shannon entropy Hq- Setting up a constrained 
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optimisation problem we derived two equilibrium equations and studied their physical implications. Ea. (l5aj) 
enforces the stability of the distribution with respect to the dynamics generated by the Hamiltonian while 
Ea. (l3bl) fixes the functional form of the distribution. Integrating the second equation we obtained a linear 
relation between Shannon entropy and the mean value of the energy, which is a strong indication that Hq 
at equilibrium has some thermodynamic properties. The physical meaning of the two equilibrium equations 
provides evidence that the maximisation of Shannon entropy is able to address the emergence of thermal 
behaviour. We also studied the relation of the proposed notion of equilibrium with quantum statistical 
mechanics. We showed that it is possible to interpret the ordinary notion of thermal equilibrium through 
the maximisation of Shannon entropy. The request to maximise the lowest among all the possible Shannon 
entropies lead to Gibbs ensembles and therefore to the quantum statistical mechanics characterisation of 
thermal equilibrium. 

Together, the physical meaning of the equilibrium equations and the proven relation between maximisa¬ 
tion of Shannon entropy and Gibbs ensemble, provide strong evidence that the proposed notion of thermal 
observable is physically relevant for the purpose of investigating the concept of thermal equilibrium. This is 
our first result. 

The result is the relevance of the proposed notion of equilibrium to address the emergence of thermal 
behaviour in a closed quantum system and especially its relation with the ETH. Using maximisation of 
Ho we were able to find a large class of observables which always thermalise and provide an algorithm to 
explicitly construct them. We call them Hamiltonian Unbiased Observables (HUOs). Using this result, we 
analysed the matrix elements of an HUO in the Hamiltonian eigenbasis to understand if they satisfy the 
ETH. The study of the diagonal matrix elements reveals that they always satisfy ETH, being constant in 
such basis. This first result on ETH has been used to prove thermalisation of the average value of HUOs, 
in connection with Reimann’s theorem about equilibration of observables. Within the picture suggested by 
Reimann, the average value of observables can equilibrate around the diagonal ensemble, if certain conditions 
hold. Arguing that these conditions are satisfied by a HUO, we applied Reimann’s theorem and obtained 
thermalisation of the average value of a HUO. 

The study of the off-diagonal matrix elements of a HUO revealed that their value is exponentially sup¬ 
pressed, in the dimension of the Hilbert space. This complete the proof that HUOs always satisfy the ETH. 
The relevance of this result for the pure-states statical mechanics is related to the two main objections usu¬ 
ally raised against ETH: the lack of predictive power for what concern both which observables satisfy ETH 
and how long it should take them to reach thermal equilibrium. The proposed notion of thermal equilibrium 
is therefore revealing its predictive power since it gives us a way of finding observables which always satisfy 
ETH, in a closed quantum system. 

We would like to conclude by putting this set of results about ETH in a more general perspective. ETH is 
one of the main paradigms to justify the applicability of statistical mechanics to closed many-body quantum 
systems. However, it is just a working hypothesis, it is not derived from a conceptually clear theoretical 
framework. For such a reason, one of the major open issues is its lack of predictability. Despite that, there 
has been a huge effort to investigate whether the ETH can be invoked to explain thermalisation in concrete 
Hamiltonian models and it’s use it is nowadays ubiqui¬ 

tous. Therefore, we think it is very important to put the ETH under a conceptually clear framework. In this 
sense, the relevance of our work resides in the fact that we obtain the ETH ansatz as a prediction, by using 
a natural starting point, the maximisation of Shannon entropy. Furthermore, using the proposed notion of 
thermal equilibrium is already giving concrete benefits. First, we now have a way of computing observables 
which satisfy the ETH and this prediction can be tested both numerically and experimentally; second, by 
studying the conditions under which local observables satisfy our equilibrium equations we are able to give 
two conditions which are necessary to guarantee the validity of ETH for all local observables. We believe 
that our investigation strongly suggest that maximisation of Hq is able to grasp the main intuition behind 
“thermalisation according to ETH” and it can be the physical principle behind the appearance of ETH. 

Further investigation in this direction is certainly needed, but we would like to conclude by suggesting a 
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way in which this new tool can be used to address the long-standing issue of the thermalisation times. From 
our investigation one can infer that Shannon entropy is a good figure of merit to study the dynamical onset 
of thermalisation in a closed quantum system. Within this picture, the time-scale at which thermalisation 
should occur for O is therefore given by the time-scale at which Ho reaches its maximum value. A prediction 
about the time-scale at which Hq is maximised will translate straightforwardly in a prediction about the 
thermalisation time. 


Methods 

Here we present the full derivation of the equilibrium equations, using the Lagrange multipliers’ technique. 
Our purpose is to find the distribution which maximises Shannon entropy of an observable O. The space on 
which such an optimisation is formulated is the space of the finite-dimensional density matrices, which is the 
space of positive semi-definite and self-adjoint matrices. The maximisation is constrained by two equations: 
the first one accounts for the normalisation of the state while the second one accounts for a fixed value of 
the average energy. 

If one wants to be absolutely rigorous, there is another constraint which needs to be imposed, which 
is the non-negativity of the density matrix. This is commonly indicated in the following way: p > 0. In 
other words, using the above-mentioned constraints there is no guarantee that, if a solution exists, it is a 
positive density matrix. This is not an actual problem, from the physical point of view, because one simply 
disregards a solution if it does not give a non-negative matrix. In such a case, the mathematical problem 
might have a solution, but not the associated physical question, which is the one we are really interested 
in. We would like to suggest an alternative route in such direction. If one wants to be completely rigorous, 
the Karush-Kuhn-Tucker (KKT) technique can be used to implement the non-negativity constraint on the 
density matrix. However, we do not think we would gain any physical insight from the use of such technique 
and therefore we are not going to explore such a route. 

We present here the derivation, in the general case of a mixed state p = 9" IV'™) (V'nl and of a 

degenerate observable O = ^ \j |j, s) {j, s|, in which {|j, s)} is a complete basis of the Hilbert space. Here 

are the two constraints: 


= Tr (p) - 1 = ^ _ 1 (27a) 

n;j,s 

CE = TripT)-Eo= Y. <lnDfsTjs,rs’D^pl-E^ (27b) 

n;j,s;j',s' 

in which 7?]”^ = {j, s|V’n)- Moreover Tjs-p's' = (j, s| T \j', s'). Exploiting Lagrange’s multipliers technique 
one dehnes an auxiliary function Aq, specific for the O observable, that can be freely optimised 


Ao[p, ^E, Aat] = Ho[p] + \nCn + XeCe 


(28) 


The derivatives with respect to the Lagrange’s multipliers Xe and Aat enforce the validity of the con¬ 
straints 


6Ao 

6Xn 

6Ao 

SXe 


= 0 Cn = 0 

(29a) 

o 

II 

tr 

o 

II 

(29b) 


while the derivatives with respect to the field variables and with respect to the statistical coefficients 
gives three equations, of which only two are independent: 
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(30) 


6Ao 


= 0 


SD>;J 


where 


SAo 

jn 

SAo 

SD jji 


-qnD):> 


-qnD):> 


log i'^qmlD 

\m,p 

log (^qm\D 


m),2 

JP I 


m),2 

JP I 


+ (1 — Xn) 
+ (1 — Xn) 


\m,p 


XeQu ,ip^. 

i,p 

XeQuY^D^^T, 


(n) 

ip 


Instead of using these equations, we use the two following independent linear combinations 


(31a) 

(31b) 


p(») SAq 


- D) 


I SAo 
SD ' 


JS 


= 0 


1 /„(n) SAq 

2 


n(») 

(/* 



= 0 


which give the two equilibrium equations that we used in the main text 


(32a) 

(32b) 


S^n{j,s) = Sn{j,s) (33a) 

- |l?g^|"log = (1 - - ^E£nU,s) (33b) 

\m,q / 
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